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ABSTRACT 


We study the early evolution of the electron fraction (or, alternatively, the 
neutron-to-proton ratio) in the region above the hot proto-neutron star formed 
after a supernova explosion. We study the way in which the electron fraction 
in this environment is set by a competition between lepton (electron, positron, 
neutrino, and antineutrino) captnre processes on free nentrons and protons and 
nnclei. Onr calculations take explicit acconnt of the effect of nnclear composition 
changes, snch as formation of alpha particles (the “alpha effect”) and the 
shifting of nnclear abnndances in nnclear statistical equilibrinm associated with 
cooling in near-adiabatic ontflow. We take detailed acconnt of the process 
of weak interaction freeze-ont in conjnnction with these nnclear composition 
changes. Onr detailed treatment shows that the alpha effect can cause signihcant 
increases in the electron fraction, while nentrino and antinentrino captnre on 
heavy nnclei tends to have a buffering effect on this qnantity. We also examine 
the effect on weak rates and the electron fraction of fluctuations in time in the 
nentrino and antineutrino energy spectra arising from hydrodynamic waves. 
Onr analysis is gnided by the Mayle & Wilson snpernova code nnmerical results 
for the neutrino energy spectra and density and velocity profiles. 


Subject headings: elementary particles - nnclear reactions, nncleosynthesis, 
abnndances - snpernovae: general 
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1. Introduction 

In this paper we examine the early evolntion of the electron fraction, in the 
post-core-bonnce supernova environment. The electron fraction is dehned to be the net 
number of electrons per baryon: 

Ye = [Ue--ne+)/ni, = l/{l + Nn/Np)-, (1) 

where rig-, ne+, and are the proper number densities of electrons, positrons, and baryons, 
respectively. The latter expression in equation (P follows from the condition of overall 
charge neutrality. Here, Nn and Np are the total proper neutron and proton densities, 
respectively, so that rii, = Nn + Np and Nn/Np is the net neutron-to-proton ratio. 

The electron fraction is a crucial determinant of nucleosynthesis produced from 
neutrino-heated ejecta in models of post core-collapse supernovae (cf., Woosley & Baron 
1992; Meyer et ah 1992; Woosley et ah 1994; Woosley & Hoffman 1992; Qian et al. 1993; 
Qian & Woosley 1996; Witti, Janka & Takahashi 1993; Hoffman et al. 1996a; Hoffman 
et al. 1996b). We follow the terminology of Fuller, Mayle, Meyer, & Wilson (1992), 
Fuller (1993) and Qian & Fuller (1995) and divide the post core bounce evolution of the 
outflowing material above the nascent neutron star into two epochs: (1) the shock reheating 
or “p-Process” epoch at times post-core-bounce tph ^ 1 s; and (2) the neutrino-driven wind 
or “r-Process” epoch at tph 1 s. 

We expect the shock reheating epoch to be characterized by chaotic outflow (Burrows, 
Hayes, & Fryxell, 1995; Miller, Wilson & Mayle 1993; Janka & Muller 1995; Herant, 
Benz, Hix, Fryer & Colgate 1994; Janka & Muller(1996)) and rapid heating by neutrino 
interactions of material behind the shock. Neutrino-heated ejecta originating from this 
epoch have been suggested to give rise to the neutron number N=50 peak r-Process material 
(Woosley Sz Hoffman 1992; Meyer et al. 1992; Woosley et al. 1994) and possibly at least 
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some of the light p-Process nuclei such as ®^Mo and ®®Ru (Fuller & Meyer 1995; Hoffman 
et ah 1996a). However, in the Woosley et ah (1994) calculations (based on the Mayle & 
Wilson supernova results) the N=50 r-Process nuclides originating in this epoch are grossly 
overproduced relative to solar system abundances. Two hxes have been proposed for the 
“N=50 overproduction problem”: (1) Fuller & Meyer (1995) invoke a modihcation of linear 
rapid outflow, a high neutrino fluence, and the alpha effect to increase W and thereby reduce 
TV = 50 overproduction; and (2) Hoffman et al. (1996a) show that as long as the electron 
fraction at this epoch can be engineered to be He ^ 0.484, N=50 overproduction is avoided. 
As a bonus, both of these hxes concomitantly suggest that some of the light p-nuclei will be 
synthesized. Hoffman et al. (1996a) hud that the light p-nuclei are produced in the correct 
proportions so long as 0.484 ^ He ^ 0.488. This epoch is characterized by relatively low 
entropy per baryon, s/k ^ 40, and relatively higher He compared to the conditions which 
may obtain at tpb ^ Is. Note that the Hoffman et al. (1996a) work implies that we may 
have to compute He to of order 1% accuracy to predict conhdently the nucleosynthesis. It 
could be that this will not be necessary, as the hydrodynamic outflow is phased in just the 
right way that a given mass element experiences the necessary He regime at the necessary 
temperature. Only future computations can address this issue. As we will show, for a 
given outflow model, predicting He histories to 1% accuracy may be next to impossible 
at this stage given our crude understanding of neutrino transport and multidimensional 
hydrodynamic effects (cf. Herant, Benz, Hix, Fryer & Colgate 1994; Janka & Miiller 1996). 

By contrast, the later r-Process epoch is characterized by considerably higher entropy, 
s/k ^ 100 — 500 (see Qian & Woosley 1996 and Meyer et al. 1992), and possibly by a well 
ordered, near steady state, outflow resembling a neutrino-driven wind (Duncan, Shapiro, 

& Wasserman 1986, Meyer et al. 1992, Qian & Woosley 1996). In fact, Woosley et 
al. (1994) have shown that the bulk of the solar systems’ r-Process material with nuclear 
mass A ^ 100 could be synthesized in this epoch. However, considerable controversy 
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surrounds the theoretical modeling of conditions in the “hot bubble” which forms in this 
regime. Though the Woosley et ah (1994) calculations yield a near perfect solar abundance 
distribution for the heavier r-Process nuclides, they are based on the very high entropy 
{s/k ~ 400) conditions obtained by the Wilson and Mayle supernova code. Not only have 
such high entropies been challenged (Qian & Woosley 1996 hnd s/k ^ 200), but even if 
the entropy were s/k ^ 300, neutrino neutral current spallation of alpha particles (Meyer 
1995) has been shown to result in a drastic and fatal reduction in the neutron-to-seed 
ratio required for the r-Process. Though models show the material in the hot bubble 
to be quite neutron rich, W ~ 0.4, Hoffman et ah (1996b) and Meyer, Brown & Luo 
(1996) have demonstrated that far lower values of T), are necessary to obtain the requisite 
neutron-to-seed ratio for the r-Process if the entropy is s/k ^ 200. At entropies this low, 
the bad effects of neutrino spallation of alpha particles (Meyer 1995) would be evaded. 
Hoffman et al. (1996b) discuss the neutron-to-seed ratio and W issues related to this epoch, 
while Fuller, Qian, & Wilson (1996) and Caldwell, Fuller & Qian (1996) discuss neutrino 
flavor mixing schemes which could give lower W and hence help the r-Process. 

In this paper we perform a complimentary study of the evolution of the electron fraction 
which concentrates on the effects of nuclear composition changes. We focus in particular on 
the early time, “low” entropy environment of the shock reheating epoch. In what follows we 
concentrate on the weak interaction balance essentially in a single outflowing mass element 
(i. e., one-dimensional outflow). We employ outflow results from the Mayle & Wilson 
supernova code to illustrate various effects bearing on W- It should be kept in mind that we 
are not predicting Y^, as the Mayle and Wilson results may not be representative of the true 
picture of supernova evolution. For example, multidimensional hydrodynamic effects could 
effectively cause different mass elements to have different time/thermodynamic histories. 
Nevertheless, we choose the simplest case (ID outflow) to elucidate the physics. 
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In section we present an overview of all of the variables which affect the calculation 
of the electron fraction. We show explicitly how the various charged current lepton capture 
processes are important. In section we discuss salient aspects of the electron neutrino 
and antineutrino capture rates on free nucleons. We explore the difference between using 
the Mayle and Wilson transport calculation-derived neutrino energy spectra as an example 
and approximate blackbody spectra. We also discuss the effect on the electron fraction 
of hydrodynamic wave-induced fluctuations in the neutrino energy spectra. In section 
we assess the role of electron and positron capture processes on W- In section we 
examine the “alpha effect” or the tendency of the formation of alpha particles to raise the 
electron fraction. We discuss the equilibrium and nonequilibrium nature of weak interaction 
freeze-out in section with particular attention to the role of nuclear composition changes 
and the role of neutrino capture on heavy nuclei. We give conclusions in section |^. 


2. Overview 

The electron fraction, which is dehned in equation ([^), can be written in the following 

way, 

w = = (2) 

i i 

Here we assume overall plasma charge neutrality, so that the sum in equation (^) runs over 
all nuclear species i with charge Zi, nuclear mass number Hj, mass fraction Xj, and number 
abundance relative to baryons Yi = Xi/Ai. Typically, the material is very hot near the 
surface of the neutron star, where essentially all of the baryons are in free nucleons. As the 
material flows away from the neutron star, it cools and alpha particles begin to form. As it 
flows farther out and cools further, heavier nuclei near the iron peak begin to form. With 
this rough evolution of abundances with radius and time in mind, we can rewrite equation 
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(0) as, 


Fe = (3) 

h 

where Xp is the mass fraction of free protons, and X^ is the mass fraction of alpha particles, 
and the summation runs over all nuclear species h heavier than alpha particles. In the 
conditions common in neutrino-heated outflow, “free” (not inside nuclei) neutrons and 
protons, alpha particles, and a few iron peak nuclei typically account for most of the 
baryons. 


The charged current weak interactions alter the electron fraction by converting 
neutrons into protons and vice versa. Most important in the region above the neutron star 
are neutrino and antineutrino capture on free nucleons and the associated reverse processes: 


z/e -t- n ^ p -f e ; (4) 

z/e + p n -1- e"*". (5) 


However, the processes of electron neutrino and antineutrino capture on heavy nuclei can 
sometimes be important in determining the overall neutron-to-proton balance (Fuller & 
Meyer 1995; McLaughlin & Fuller 1995): 


z/e + A(Z,N)- 

A(Z -f-1, N — 1) -|- e ; 

(6) 

Pe + A(Z,N)- 

A(Z-l,N + l) + e+. 

( 7 ) 


In these expressions. A, Z, and N are the total number of nucleons, proton number and 
neutron number of the nucleus, respectively. The reverse rates of electron and positron 
capture on heavy nuclei are generally negligible for the conditions in which these nuclei 
form. The ratio of neutrons-to-protons and Yg in neutrino-heated material flowing away 
from the neutron star is set by a competition between the rates of the processes in equations 
(I) and (^ and (|^) and (^ and the overall material expansion rate (or outflow rate). In 


fact, it has been shown that where the rates of these processes are rapid compared to the 
ontflow rate, a characteristic weak steady state or “chemical eqnilibrinm” obtains (Qian et 
ah 1993; Qian 1993). The weak freeze-out radius is dehned to be where the rate of n ^ p 
interconversion as set by the rates of the processes in equations i). (I). i) and d^) falls 
below the material outflow rate. Though the details are complicated, Y^, at small radius is 
set principally by the processes in equations (^) and (j^); whereas, at larger radius and at 
later times, the processes in and (|^) can also make a contribution. 

The rate of change of the electron fraction of an outflowing mass element resulting 
from weak interactions may be written as follows (this is a generalization of the treatment 
in Qian et ah 1993): 

~ + ^e+)Xn — (Ai/g -|- Ae-)Xp-|-^ ~ ~ ^he-), (8) 

where v{r) is the radial velocity field above the neutron star, and f is a time development 
parameter. In equation (^, the sum on h indicates a sum over all heavy nuclei. Here 
Xn is the mass fraction of free neutrons, while Xh is the appropriate mass fraction of 
a heavy nucleus with mass number A^. We dehne and Xp^ to be the forward rates 
of the processes in equations (^ and (^, respectively. Similarly, X^,- and Ae+ are the 
reverse rates of these processes, respectively. Finally, XhuX^he+) and XhvX^he-) represent 
the electron neutrino (positron) and antineutrino (electron) capture rates, respectively, on 
the nucleus with index h. Note that equation (8) contains all the effects which bear on 
Ye and the neutron-to-proton ratio in a given mass element. Hydrodynamic motion can 
influence this weak balance via the position of the weak freeze-out radius (Qian et al. 1993) 
which is determined largely by the velocity held v{r) (more properly, in multidimensional 
hydrodynamic regimes we should write v{r, 6, 4>,t))- 

Near the neutron star, before the iron peak nuclei form, the last term in equation (8) 
is zero. The neutrino and antineutrino capture reactions on free nucleons are usually the 
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fastest of the nuclear charged current weak interaction processes. The exception occurs 
in the region very close to the neutron star, where electron and positron capture on free 
nucleons become comparable to the corresponding neutrino capture rates. In the region 
near the neutron star where the electrons are relativistically degenerate, Ag- can be large 


and can essentially set to quite low values. By contrast, at late times and/or large radius. 


the free nucleons become absorbed into nuclei, and the last term in equation (|) becomes 
larger than or competitive with the other two terms (in this environment, Xhe+ and Xhe- are 
always small compared to Xhp^ and/or Xh^J- Here we include no term which is proportional 
to Xa, since charged current capture rates on alpha particles are too small to influence the 
electron fraction in the relevant conditions. 

We can give a formal solution to equation (|]) in the limit where there are free nucleons 
and alpha particles only, and no heavy nuclei: 



+ [-^q/2 + Xn /(Ap + An)] J 



In this expression the initial conditions are denoted by the subscript i and the hnal 
conditions are denoted by the subscript /. Here, An = Ae+ + A^^ and Xp = Ag- + Ap^ are the 
total neutron and proton destruction rates, respectively, resulting from weak interaction 
processes. When the system is in weak equilibrium, the rate of change of the conditions 
with time or radius in the plasma will be slow compared with the magnitude of the weak 
rates A„ and Xp. The hrst term in equation (j^) will be negligible when weak equilibrium 
obtains. The second term in this equation will be the value of the electron fraction as 
the weak equilibrium limit is approached. The exponential factor in the integrand of the 
last term is very small except for t near tf. This factor will not be small for integration 
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intervals, t — tf = 6t l/(An + Xp). As long as the product of the time derivatives of 
and Xn/{Xp + A„) with 6t is small, the last term in equation (^) will also be small. This 
condition is usually satished where weak equilibrium obtains and where the alpha particle 
mass fraction is only slowly varying with time. If the system is not in weak equilibrium, 
and/or the alpha particle mass fraction is changing signihcantly with time, then in principle, 
all of the terms in equation (j^) may be necessary for calculating the electron fraction. 

If weak equilibrium is established, alpha particles are absent, and the electron and 
positron capture rates, Ae- and Ae+, are both zero, then the electron fraction will be, 

y; = 1/(1 + KjK.) = y;„. (lo) 

This will be a reasonable approximation to the true electron faction at about tp/, ~ 1 s in the 
region just below where alpha particles form (Qian et al. 1993). We shall use the quantity 
Wo as a hrst estimate of the electron fraction. We discuss computational estimates of Wo in 
detail in section |^. 

Although Wo provides a good hrst estimate of the electron fraction, the electron 
and positron capture rates on free nucleons can also have some inhuence. In fact, at a 
plasma temperature of ~ 2 MeV, the electron and positron capture rates can make a 
signihcant contribution in the determination of the electron fraction. If the system is in 
weak equilibrium and alpha particles are absent, then the electron fraction in this regime of 
high plasma temperature is more accurately given by, 

W = 1/[1 + {Xpe + Ae-)/(Ae+ + Ajy^)] = (1 + Xp/Xn) (11) 

This can be seen from equation (^ above. At such high plasma temperature, ~ 0 
usually will be a good approximation for the relevant entropies in the supernova models we 
consider. The effect of electron and positron capture on W will be discussed in section ^ 
We generalize this discussion to include the Xa 7 ^ 0 case in section 
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Fuller & Meyer (1995) first pointed out that the mass fraction of alpha particles can 
influence Y^. This effect can be readily gleaned from equation (j^). When the alpha mass 
fraction rises (plasma temperature ~ 0.5 MeV), free nucleons are typically absorbed into 
alpha particles. Each alpha particle removes an equal number of free protons (2) and free 
neutrons (2) when it forms. If the outflowing material in the supernova is neutron-rich 
before the alpha particle formation, then as the alpha particle mass fraction increases the 
residual free nucleon gas will tend to be enriched in neutrons. Since the charged current 
capture rates on alpha particles are negligible in this situation, the only interactions 
that change the electron fraction are lepton captures on free nucleons. However, since 
alpha particle formation has left mostly free neutrons, the neutrino and positron capture 
interactions will have the effect of turning some of these “leftover” free neutrons into 
protons. In turn, this will cause the total electron fraction to rise. We follow Fuller & Meyer 
(1995) and call this increase in W due to alpha particle formation and weak interactions the 
“alpha effect.” We give a detailed treatment of this effect in section H. 

Since the material may not be in weak equilibrium at the time of nuclear reaction rate 
freeze-out and nucleosynthesis, it is necessary to examine the effect on the electron fraction 
of slow freeze out from weak equilibrium. This is done in section |^, where we present the 
results of a nonequilibrium calculation of W- The influence on W of neutrino captures on 
heavy nuclei will also be discussed in section As is evident from equation (H), these 
processes in principle can impact the hnal value of the electron fraction. Although it will 
turn out that the material is no longer in weak equilibrium at the time of formation of 
iron peak nuclei in the models we consider, the neutrino captures on heavy nuclei (plasma 
temperature ~ 0.3 MeV) still have some influence on W- 
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3. Neutrino and Antineutrino Capture Rates on Free Nucleons 


In this section we compare the relative magnitnde and effects on Y^, of and Ap^. For 
the calculation of these rates we follow the prescription of Fuller & Meyer (1995; 1996) and 
McLaughlin & Fuller (1995; 1996). In computing these rates, we hrst assume a blackbody 
distribution for the neutrino spectrum, which we normalize by the neutrino luminosity (see 
Qian & Fuller 1995). In Figure we depict a zero chemical potential blackbody neutrino 
(z/g) energy spectrum (curve without circles) with average neutrino energy and normalizing 
luminosity taken from the numerical results of the Mayle and Wilson supernova calculations 
at tpb ~ 0.575 s“^. Also shown in this hgure is the actual neutrino energy spectrum (circles) 
at large radius at tpf, ^ 0.575 s“^ as given by the detailed transport calculations of Mayle 
and Wilson. Of course, the actual neutrino energy spectra may differ considerably between 
different models which employ different transport calculations. We here use the Mayle and 
Wilson results only as an example. 


In our calculations of and Ap^ we generally assume that the material is far from 
the neutron star, so that the distance dependence of the rates is simply oc where rj 
is the distance from the center of the neutron star in units of lO^cm. We further assume 
that the hnal state electrons are very relativistic, so that the rate may be expressed with 
Fermi integrals. In evaluating the Fermi intergals for the antineutrino capture rates, we use 
the approximation {rup — rrin — me)/Tp -C 0, where is the mass of the neutron, nip is 
the mass of the proton, rUg is the mass of the electron, and Tp is the temperature of the 
antineutrino distribution function (we assume zero chemical potential). The capture rates 
are then approximately. 


A;,^ ^ (0.1945 s 


-1 


( K \ (T^ 

yi0®^ergs“^y VMeVy [vj 


4|Ci; 


( 12 ) 


Ap^ (0.2000 s ) exp 


-1.804 MeV 


Lp 


T- 


) (l0®^ergs“^y VMeVy 


Co. 


(13) 
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In this expression we have approximated the initial and hnal state lepton kinematics as 
completely relativistic, and we have neglected hnal state lepton blocking (McLaughlin & 
Fuller 1995). The terms Ci and C 2 are dehned in the following way: 

Cl ^ 1 + (0.6283 MeV)/T^ + (0.1292 (14) 

C 2 -I + (1.158 MeV)/T^ + (0.600 MeY‘^)/T^ + (0.1778 MeY^)/Tl (15) 

These estimates demonstrate the dependence of the rates on neutrino (antineutrino) 
luminosity, the temperature, and distance from the neutron star center, rj. 

The rates Ap^ and are approximately proportional to the antineutrino and neutrino 
temperature, respectively, when these temperatures are large. However, at lower neutrino 
and antineutrino temperature they have a more complicated temperature dependence 
through the terms Ci and C 2 and the exponential term in equation ([T^). The exponential 
term in the antineutrino capture rate has its origin in the energy threshold, nin — mp + me- 
Of course, there is not a corresponding term in the neutrino capture rate, since neutrino 
capture on free neutrons has no energy threshold. We estimate the temperature of the 
blackbody distribution to be ~ E/3.15, where E is the average energy of the neutrino 
or antineutrino spectrum as appropriate. At tpb ~ 0.500 s, the neutrino temperature is 
~ 3.85 MeV and the antineutrino temperature is Tp 4.93 MeV in the Mayle and Wilson 
numerical calculations. 

The ratio of the antineutrino and neutrino capture rates on free nucleons is 

^ ^ 1.029 (p) exp [(-1.804 MeV)/Tp]Ji, (16) 

where the term Ji is dehned as, 

Ji = ^ 1 + (1.158 MeV)/Tp - (0.6283 MeV)/T^. (17) 

Since the neutrino and antineutrino temperatures are relatively high, we have ignored terms 
of higher order in these temperatures in Ji. Because of the neutron excess in the neutron 
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star, Tp^ will be generally larger than everywhere above the neutrino sphere (we ignore 
possible neutrino flavor mixing effects - see Qian et ah 1993, and Qian & Fuller 1995). Also, 
at this epoch, the total energy luminosity of the neutrinos is about 10% smaller than that of 
the antineutrinos in the Mayle and Wilson supernova results. In order for the antineutrino 
capture rate to dominate over the neutrino capture rate, the antineutrino temperature 
and luminosity must be sufficiently greater than the corresponding neutrino temperature 
and luminosity to overcome the effect of the energy threshold. Since the neutrino and 
antineutrino capture rates have the same dependence on distance from the neutron star, the 
relative importance of the rates does not change explicitly with radius in this formulation. 

The z/e and Ug temperatures and luminosities change with time, and therefore the 
neutrino and antineutrino capture rates are changing implicitly as the material flows away 
from the neutron star. The ratio of our approximate antineutrino capture rate on protons 
to our approximate neutrino capture rate on neutrons has been plotted as a function of time 
(or epoch) in Figure |^. The antineutrino capture rate is always larger that the neutrino 
capture rate in the indicated time interval, and this difference is increasing with time. The 
curve in Figure § has been smoothed into a straight line for illustrative purposes. Variations 
of this ratio from a strictly increasing function will be discussed further below. 

Since the ratio Xp^/Xu^ is important for the calculation of the electron fraction, it is 
better to use the numerically computed energy spectrum to calculate the neutrino and 
antineutrino capture rates. An example of such a numerically calculated energy spectrum 
taken from the results of the Wilson and Mayle code is shown in Figure The circles 
in Figure show the points produced by the detailed neutrino transport computations of 
this code. In order to calculate the capture rates from these numerical points it is first 
necessary to interpolate to And a value of the energy occupation probability for all neutrino 
energies. Given such an interpolation scheme, the capture rates can be calculated with the 
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expressions, 



(18) 


Here, is the interpolated valne for the nnniber flnx of (anti)nentrinos at energy E. 

The (anti)nentrino sphere radins is given by Ru{p), and r is the distance from the nentron 
star center at which the rate is to be calcnlated. The distance from the nentron star center 
at which the nentrino spectrnm has been evalnated is tq. In the case of the Wilson and 
Mayle code, this distance is ro = 3 x lO^cm. Assnming relativistic lepton kinematics, and 
no hnal state lepton blocking, the cross section at energy E is given by 



(19) 


where Q is the approximate nnclear Q-valne (see eqnation 4 in McLanghlin & Fnller 1995). 
The Q-valne for nentrino captnre on nentrons is Q ~ 1.293 MeV, while for antinentrino 
captnre on protons it is Q ~ —1.293 MeV. The energy threshold is Eth = 0 for nentrino 
captnre on nentrons while it is Eth = (^^n — + ’^e) ~ 1.804 MeV for antinentrino 

captnre on protons. 

Table |I] illnstrates the discrepancies in captnre rates and electron fractions calcnlated by 
different methods. In the hrst colnmn of this table, the method of calcnlating the nentrino 
or antinentrino spectrnm is given. In the second and third colnmns, the corresponding 
nentrino and antinentrino captnre rates on free nentrons and protons, respectively, are 
given. In the fonrth colnmn, the qnantity, Wo = 1/(1 + is tabnlated. All rates are 

calcnlated at a distance from the nentron star center of r = 4.793 x 10® cm. All qnantities 
are calcnlated at time tpb = 0.575 s in the Mayle and Wilson calcnlations. The 


temperatnres employed for the blackbody spectra in the hrst row are = 3.724 MeV and 
Tp^ = 4.835 MeV. 


The hrst row of Table shows the resnlts of the calcnlations with zero chemical 
potential blackbody nentrino and antinentrino energy spectra. Snch a blackbody energy 
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distribution tends to overestimate the (anti)neutrino capture rate. This is due to depletion 
in the high energy tail of the actual transport calculation-derived neutrino energy spectrum 
relative to the zero chemical potential blackbody approximation. This depletion is caused 
by (anti)neutrino absorption and scattering above the neutrino sphere. The blackbody 
spectrum tends to underestimate the electron fraction, as shown in the column labeled Yeg. 


The second row of Table |I| gives the resulting capture rates and electron fraction when 
cubic spline interpolation is used to generate the segments of curve between the points 


on the numerically-derived (anti)neutrino energy spectra, and equations (p!8[) and (|T^ are 
used to calculate the capture rates. In order to check the suitability of this method for our 
purposes, we tried another interpolation scheme to estimate the capture rates. We £t curves 
between the numerically-derived spectrum points with pieces of blackbody distribution. 
The Mayle and Wilson spectrum points themselves are derived from the numerical transport 
calculations and are produced with an interpolation procedure in some ways similar to that 
employed here. 


Each such blackbody segment had a different temperature. As the neutrinos leave the 
neutron star, different energy neutrinos decouple at different radii and temperature. This 
motivated our htting of the neutrino spectrum by pieces of different thermal distributions. 
The results are shown in the third row of Table ||. The difference between the approximate 
electron fraction Wo obtained in this last case and in the case of cubic spline interpolation 
was about 0.3%. The difference in Wo between the case with cubic spline interpolation and 
the case with a single zero chemical potential blackbody distribution was about 4%. Since 
it may be necessary to obtain the value of the electron fraction to within 1% (Hoffman et 
ah 1996a; Qian & Woosley 1996), cubic spline interpolation gives an adequate result. This 
is the method that we have employed in our calculations of the neutrino and antineutrino 
capture rates throughout the remainder of this paper. 
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Figure ^ shows the value of the parameter at different times, calculated using 
cubic spline interpolation for the neutrino and antineutrino energy spectra taken from 
a particular run of the Wilson & Mayle supernova code (crosses). There is signihcant 
variation in this quantity during the time increment {At ~ 0.025 s“^) between typical 
snapshots of the neutrino and antineutrino energy spectra at tph ^ 0.575 s. In the hot 
bubble above the neutrino sphere, hydrodynamical waves are propagating back and forth. 
Spherical convergence may amplify these waves near the neutrino sphere. Multi-dimensional 
non-spherical geometries may complicate this picture considerably (Burrows et ah 1995; 
■Ta nk a & Muller 1995). The variations evident in Figure ^ may be due to these physical 
effects. The value of Wo obtained with a similar procedure, but employing a different 
Mayle and Wilson numerical model for this epoch, suggests smaller variations. In any case, 
since we may need to compute W to of order ~ 1% accuracy for nucleosynthesis purposes 
(Hoffman et al. 1996a), it is important to be aware of this issue. Since this variation in 
Wo is relatively large ~ 3%, it will be necessary to resolve this and other issues before the 
electron fraction can be reliably calculated to accuracies of approximately ~ 1% during 
this epoch. It is obvious from these considerations that other supernova outflow models, 
employing different (anti)neutrino energy spectra and/or different hydrodynamic conditions 
may well give a W very different from that computed here. Our results are meant to be 
only illustrative of the physics required to estimate W to 1% for any model. 

For illustrative purposes we have chosen to draw a straight line through the time period 
0.575 s — 0.700 s (during which the point-to-point variations in Wo small) yielding the 
linear relation. 

Wo ~ 0.4606 - 0.7120(t - 0.5750). (20) 

This line, which is shown in Figure is only meant to demonstrate the effect of the 
increasing disparity between the neutrino and antineutrino energy spectra with time on Wo- 
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The neutrino and antineutrino energy spectral data that we use here conies from 
snapshots taken from the Mayle & Wilson output at various times post-bounce. At each 
time, these spectra have has been evaluated at a radius of 3 x lO^cm. However, the neutrino 
and antineutrino energy spectra will change at different points in the region above the 
neutron star due to weak interaction processes in the plasma. For the situation we are 
considering, we are far from the neutrino sphere, and close to the radius tq of spectral 
quantity evaluation. However, it is clear that the effects of neutrino emission, absorption, 
and scattering must also be included for a precise determination of Wo in fh® region near 
the neutron star. We believe that these effects are smaller than the observed oscillation in 
the quantity Wo apparently caused by local variations from hydrodynamic waves. 


4. Electron and Positron Capture Rates on Free Nucleons 

The processes of electron and positron capture on free nucleons are given in equations 
(^) and (|^). These processes are important in regions close to the neutron star. They play 
a role in determining the position of the gain radius. Neutrino interactions cause a net 
positive heating of the region above the gain radius (Bethe & Wilson 1985). At an early 
epoch (tpb ^ Is), the weak freeze out radius, or distance from the neutron star center at 
which the characteristic rate of weak interactions becomes smaller than the rate of material 
outflow in the supernova, occurs sufficiently beyond the gain radius so that electron and 
positron capture rates on free nucleons have fallen considerably below the corresponding 
neutrino and antineutrino capture rates. At this distance, and at all points further away 
from the neutron star, the electron (positron) capture rates are always much smaller than 
the (anti)neutrino capture rates on free nucleons. 

Here we give estimates of the electron and positron capture rates on free nucleons, 
analogous to the estimates of neutrino and antineutrino capture rates presented in section ^ 
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(equations ^ and |^. We employ Fermi-Dirac distributions for the electrons and positrons, 
since these particles are well approximated as being in thermal equilibrium with the plasma. 
This is a much different situation than the case of neutrinos, since the neutrinos decouple 
from the plasma very close to the neutron star. In producing these estimates, we have 
assumed that the electrons are relativistic, so that the appropriate phase space factors 
may be reduced to Fermi integrals (Fuller, Fowler, & Newman 1985; 1982ab; 1980; Fuller 
1982). Although the electrons and positrons are not completely relativistic throughout the 
entire period of interest for nucleosynthesis, their capture rates on free nucleons are only 
important at fairly high temperatures (Tg > IMeV) and relatively small distances from the 
neutron star. Therefore, this approximation is sufficient to demonstrate the general effect 
of these capture rates on the electron fraction. The estimates of the rates are given by. 


Ae- ^ (1.578 X 10"2s“M exp ( 

KmeC^J V 


Ae+ ^ (1.578 X 10“^ s-^) 




rripC^ 


exp 


'-1.293 + pe- 
T 

e 

-0.511 -Pe- 



( 21 ) 

Q. 

( 22 ) 


Here C 3 and C 4 are dehned in the following way: 


(Fs ^ 1 + (0.646 MeV)/Te + (0.128MeV^)/Tf; (23) 

(^4 ^ 1 + (1.16 MeV)/Te + (0.601 MeV^)/T^ + (0.178 MeV^)/T^^ + (0.035 MeV^)/T^^ (24) 

In these expressions, as everywhere in this paper, the total electron chemical potential Pe- 
is dehned as the kinetic chemical potential plus the rest mass of the electron. We have 
estimated the value of the Fermi integrals (see Fuller, Fowler, & Newman 1985 equations 
5a-f) assuming {nip — nin + fie-)/Te « 0. Clearly, this assumption is only valid when the 
electrons are not very degenerate. The region we are studying is far from the neutron star, 
and at sufficiently high temperature and low density that this approximation is reasonable. 
These rates depend strongly on temperature, which causes their magnitude to fall quickly 
in the outhowing material. Because of the many approximations involved in obtaining 
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these expressions, they may not faithfully represent the true rates to the accuracy necessary 
for nucleosynthesis calculations. Therefore, we use them only to illustrate their influence 
on the electron fraction. For example, convection or other multidimensional effects could 
necessitate following the lepton capture rates in degenerate conditions which would modify 
the evaluation of the rates in equations (21) and (22). The non-equilibrium calculation 
presented in section ^ includes a more precise calculation of the electron and positron 
capture rates. 


The hrst question that we address concerns the relative magnitude of these rates. We 
can combine the approximate expressions for the rates into a ratio: 


^ ss exp [(-0.782 + 2p,-) /TJ (^) ; 


(25) 


where we intend Pe- and Tg to be expressed in units of MeV. The reduction in the 
ratio caused by the threshold for the electron capture reaction is contained in the term 
exp(—0.782 MeV/Tg). The expression in equation (E^ also demonstrates the effect of the 
chemical potential in reducing the positron capture rate. In general, the temperature in the 
region of interest is not high enough to neglect the additional multiplicative term, 
which comes from the evaluation of the phase space factors. In Figure ^ we show a plot of 
this ratio. The curve Ag-/Ag+ is depicted starting (far left) in conditions where the system 
is in weak equilibrium, at a temperature of about Tg 2 MeV. This corresponds to a time 
of tpb = 0.575 s in the hgure. 


In our calculated estimates of Ag-/Ag+ in Figure ^ we have utilized an initial Vg from 
the weak equilibrium condition, a time dependent density £t to the Mayle and Wilson 
supernova computational results, and a constant value of the entropy per baryon set to 
s/fc ~ 80. In an actual supernova we would not expect the entropy of the outflow to be 
constant, nor would we expect s/k = 80 to obtain necessarily. We pick s/fc ~ 80 and 
approximate the outflow as adiabatic for illustrative purposes only. In section 6, we consider 
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other values of s/k. Our time dependent density fit can be taken to represent the density 
history of an outflowing Lagrangian mass zone during a limited time period. The functional 
form we have adopted for this £t is: 

p ^ exp (-310.154 + 1811.43t - 3215.62^^ + 1839.84^^) g cm'^; (26) 


where t is time post core bounce. We have calculated all other thermodynamic variables, 
such as the temperature and chemical potential, from this density fit and assumed entropy. 
If the assumed entropy is decreased, then the effect of electron degeneracy will exhibit 
itself in the ratio of the rates. In this case, the ratio of the electron capture rate to the 
positron capture rate will be greater than that shown in Figure]^, since the number density 
of electrons will be enhanced, while the the number density of positrons will be suppressed. 
The ratio Ae-/Ae+ decreases at lower temperature (larger time), due to the increasing 
importance of the threshold for electron capture, and also because of the increasing ratio of 
positron number density to electron number density. It can be seen from this figure that the 
ratio of the electron capture rate to the positron capture rate exhibits much more variation 
than the ratio of the antineutrino capture rate to the neutrino capture rate on free nucleons. 


The ratio of the positron capture rate to the neutrino capture rate on free nucleons 
can serve to illustrate an important point about the evolution of conditions with radius and 
time in the outflowing material. The ratios may be expressed as follows: 


A. 


2.34 


( Te 2 AO^^ergs^ /^MeV\ ^ fCi\ 


where the electron rest mass is rUe ~ 0.511 MeV. In Figure we plot \e+l^v^ as a function 
of time. It is evident from this hgure that as a mass element flows away from the neutron 
star, Ae+ drops off more quickly than does The drop in \^+ is a result of the rapid 
decrease in plasma temperature with radius, while the drop in X^^ with radius simply 
reflects the relatively slower fall off in the neutrino flux (the l/r^ term in equation |^. The 
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decrease in the electron chemical potential /ig- with radius has a relatively smaller effect in 
comparison to the fall in plasma temperature in setting Xe+Z^^u^- 

Since, as shown above, the electron and positron capture rates on free nucleons are 
small relative to the corresponding neutrino and antineutrino capture rates, the electron 
fraction may be written as an expansion in the small parameter Xe+/X^^. Employing this 
small parameter, we can expand the expression for the electron fraction in equation (P) to 
yield. 



( 28 ) 


The first order term in this equation is sufficient to show the effect of electron (positron) 
capture rates on the “equilibrium” T).- Since the lepton capture processes are not necessarily 
in true chemical or steady state equilibrium, by “equilibrium” Ye here we mean the Ye that 
would obtain if the system were in true equilibrium. As the material flows outward in the 
region of interest, the term Xe+jXy^ decreases rapidly. This term represents the competition 
between increasing distance from the neutron star and decreasing plasma temperature. 
However, the term [1 — (1 + Ae-/Ae+)/(l + Aiy^/A^^)] actually shows a slight increase with 
radius or time. This is caused by the decreasing ratio \e-/Xe+- In the example discussed 
here the increase in the term [1 — (1 + Ae-/Ae+)/(l + Ai,^/ApJ] is overwhelmed by the 
decreasing positron capture rate in the leading Xe+/Xy^ term, so that the change in the 
“equilibrium” electron fraction with radius is dominated by the behavior of Xe+/Xy^. 

Figure ^ shows this “equilibrium” electron fraction. The upper curve labeled 
“equilibrium T).” includes the effect of electron and positron captures on free nucleons. The 
lower curve in this figure labeled “Hep” includes only neutrino captures on free nucleons. 
At high temperature (Tg ~ 2 MeV) corresponding to early times, the electron and positron 
captures make a considerable difference in the electron fraction (on the order of ~ 10%). 
In fact, we have slightly underestimated their effect by using the approximate expressions 
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for the rates in equations (pi]) and However, it is clear that these processes become 

less important with time, as the “equilibrium” electron fraction slowly asymptotes to the 
value of Heo- more electron degenerate conditions, the electron capture rate may be 
greater than the positron capture rate, causing the “equilibrium” electron fraction to lie 
below Heo at early times. In this case, the electron and positron capture processes change 
the “equilibrium” electron fraction by about 2% just before alpha particle formation. We 
emphasize that this conclusion depends on our particular Mayle & Wilson outflow history. 
Different models will give disparate W values (the results of section 6 for different entropies 
could be used to gauge how different Yj. might be in other, more realistic models). 


5. The Alpha Effect 


In this section we use the “equilibrium” W to discuss qualitatively the behavior of the 
electron fraction during alpha particle freeze-out. The formation of alpha particles occurs 
while the outflowing plasma is still in nuclear statistical equilibrium (NSE). Therefore, the 
mass fraction of alpha particles at a given density and temperature is well approximated by 
the Saha equation: 


X{N,Z) = 

pNa 


2TTh^NA 

k%. 


3(A-l)/2 


/pNaX^Y (pNaX, 


V A„ 




\ 




3.256 X l')^;f2;f2g,,p(28.29MeV/r,), 


exp[(3„(iV,Z)/J;rJ 

(29) 

(30) 


(Te/MeV)' 

where X{Z, N) is the mass fraction of the nuclear species with Z protons and N neutrons, 
G{Z, N) is the partition function, p is density, Na is Avagadros’ number, k is the Boltzmann 
constant, Qn{N,Z) is the nuclear Q-value, and An and Ap are the neutron and proton 
atomic masses, respectively. In equation Pio is the density in units of 10^°gcm 
In order to estimate the number density of alpha particles, we use a plasma temperature 
calculated from a given entropy and the htted density function (equation as discussed 
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in the last section. 

To zero order in the ratio Xe+Z^v^, the electron fraction in the presence of an alpha 
particle component is given by (see equation 

~ ^a(0.5 — Teg). (31) 

The magnitude of the “alpha effect” on T). is proportional to the number of alpha particles 
and also to the difference between T);g and 1/2. If there are equal numbers of neutrons and 
protons, then there is no change in the electron fraction due to alpha particle formation. 
However, if the neutron-to-proton ratio is not unity, then the effect of alpha particle 
formation will be to drive T/ closer to 1/2. The expression in equation (|3T|) does not include 
the small effect of electron and positron capture during this period. These captures may be 
accounted for by including the terms to hrst order in the ratio \e+/Xu^\ 


Ye ~ T/g 

+^eoX~ [1 ~ (1 + Xe-/Xe+)/ (1 + X^^/ AjjJ] 

+Xq( 0.5 — T/g) 

(w) fw -1 -1) ( a"a"+i ) 


(32) 


The second line in this expression is the correction to the “equilibrium” Ye from electron 
and positron captures on free nucleons discussed in the last section. The last line in 
equation (|3^) contains an additional correction term for the “equilibrium” T/ due to 
the electron and positron capture processes on free nucleons when alpha particles are 
present. These additional terms give a more complete description of the “equilibrium” T/ 
when alpha particles are present. Figure shows the value of l/g and the value of the 
“equilibrium” Ye as a function of time for our example outflow trajectory. At early time and 
high temperature, the influence of the electrons and positrons can be seen by the increase 
in the “equilibrium” Ye over T/g. The minimum in the “equilibrium” Ye curve occurs at 
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the point where alpha particles begin to form. The increase in “equilibrium” Y^, as time 
increases subsequent to the minimum is due to the “alpha effect.” Since the true Y^. is not 
in equilibrium at the time of alpha particle formation, the increase in the true Ye will be 
much smaller than the increase in the “equilibrium” Ye shown in the hgure. The actual 
increase in the true Ye depends on the rate of alpha particle formation. 

Since the process of incorporating nucleons into alpha particles is clearly an important 
factor for determining the electron fraction, we wish to explore the relationship between the 
rates of neutrino and antineutrino capture on free nucleons and the rate of alpha particle 
formation. Assuming a slow and constant rate of alpha particle formation = Ca, and 
assuming constant neutrino and antineutrino capture rates, we can integrate the differential 
equation for Ye (equation |^) to obtain: 

n « Ye, + (i/(2n,) -1). (33) 

Here Ye^ is the value of the electron fraction before alpha particle formation begins. 
Equation (^) will be valid during the epoch where the alpha particle abundance is changing 
and is accompanied by a continuously changing electron fraction. The presence of the term 
in equation (|33D shows that there is a competition between alpha particle formation 
and the neutrino capture on neutrons in setting Ye. Either a faster neutrino capture rate or 
a slower rate of alpha particle formation will result in an increased change in the electron 
fraction. If the material is neutron rich this change will be positive. The last term in this 
expression shows that there will be little effect from alpha particle formation when the 
electron fraction is close to Ye = 0.5. 

The change in the electron fraction due to alpha particle formation may be estimated 
from equation (PSI). At tph ~ 0.5s the appropriate quantities taken from our calculations 
give Xp^/Ca ~ 5s“^/50s“^ ~ 0.1 and (l/(2Feo) — 1) ~ 0.3. Therefore, the change in the 
electron fraction as computed from equation (|33D will be ~ 0.003, approximately 1%. 
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The non-equilibrium calculation presented in section gives the results of a more exact 
and detailed treatment of these issues. These results can be used to estimate how much 
different would be for different models of thermodynamic history and (anti)neutrino 
energy spectra. 


6. Non-Equilibrium Calculations 

In this section we focus on the non-equilibrium aspect of the calculation of the electron 
fraction. In the previous sections we have focused on the “equilibrium” Ye, in order to show 
the general trends caused by each of the variables which influence the electron fraction. 
However, the relative magnitudes and effects of these factors can only be determined 
dehnitively by a full non-equilibrium calculation. 

For this calculation, we employ an NSE computer code to keep track of thermodynamic 
variables, weak reaction rates and the electron fraction in a representative mass element 
of outflowing material. As input, the code utilizes the density £t from equation (^), an 
assumed initial electron fraction and a constant entropy as before. The outflow velocity of 
the material is obtained by interpolating with cubic splines between the velocity given in a 
particular run of the Mayle and Wilson results at different time slices. Other thermodynamic 
variables, such as temperature and chemical potential are calculated self-consistently from 
these parameters as outlined in section but now modified where appropriate to include 
the effects of heavy nuclei. We calculate the relative abundances of free nucleons, alpha 
particles, and heavy nuclei from the thermodynamic variables and the current value of the 
electron fraction. The electron fraction is updated at each time step in our calculation by 
following all the charged current weak reaction rates. 

During the periods when only free nucleons are present, the electron fraction determines 
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the values of Xp and X^- After alpha particles form, relative abundances are calculated 
with the electron fraction and the nuclear Saha equation for alpha particles, equation 
(17). During the periods when heavy nuclei are present, we use the liquid drop model and 
the prescription of Fuller (1982), Bethe, Brown, Applegate and Lattimer (1979; hereafter 
BBAL), Baym, Bethe, and Pethick (1971) and Lamb et ah (1978) to calculate the neutron 
and proton chemical potentials. In this model, the total energy of the nucleus with mass 
number A can be approximated as, 

Wn ^ IFbuik + + IFcouiA^A. (34) 

In this expression, BAuik is the energy of bulk nuclear matter, while hFsurf and hFcoui are the 
coefficients of the surface energy and coulomb energy terms, respectively. The values for 
these coefficients given in BBAL and employed by Fuller (1982) are, 

IFsurf ~ 290(Z/A)2(1 - Z/Af- (35) 

hFcoui ~ 0.75(Z/A)2(1 - 0.236p}^^ + 0.00194pi2), (36) 

where pi 2 is the density in units of lO^^gcm”^, and Z is the number of protons in the 
nucleus. The value of the mean nuclear mass in NSE is obtained by minimizing BAr, 
yielding the condition that the nuclear surface energy should be twice the Coulomb energy, 
from which we derive, 

A ^ 194(1 - (Z/A)f(l - 0.236pi^^)"^ (37) 

The difference in neutron and proton chemical potentials is given by, 

= f^n — l^p ^ 250 [0.5 — [ZjA)\ — WsurfA ^ 1(2’ jA) \y — [Z /A)] | ’ 

The neutron chemical potential in this scheme is (Fuller 1982; BBAL), 






16 + 125[0.5 - (Z/A)] - 125(0.5 - {ZlA)f 


IF,„rf (3 - 7(Z/A) 

AV3 2[1 - (Z/A)]‘ 


(39) 



For a recent calculation and discussion of the nuclear equation of state see Lattimer & 
Swesty (1991). 


The neutron chemical potential has the same value for the neutrons inside the nucleus 
and the neutrons in the free neutron gas. In the dilute Maxwell-Boltzman limit which 
obtains for free nucleons, the mass fraction of free neutrons will be, 

.Y„ « 79 PVMeV)y^ exp{i^JZ), (40) 

PlO 

while the corresponding mass fraction for free protons will be, 

Xp ^ X„exp(-/i/Te). (41) 


The mass fractions of heavy nuclei are calculated by using the nuclear Saha equation. 
The partition functions and binding energies employed in this prescription were taken from 
Woosley, Fowler, Holmes and Zimmerman (1978). Only heavy nuclei which contribute more 
than a few percent to the total heavy nucleus mass fraction at a given time are retained, 
and the abundance distribution is normalized to insure that 


A'p + A'„ + A'„ + J;A'j = 1. 


(42) 


We allow the electron fraction to change by calculating weak capture rates on free 
nucleons and heavy nuclei. In the case of (anti)neutrino capture on free nucleons, the 
neutrino distribution functions (smoothed with time) are employed as discussed in section 
The electron and positron capture rates are calculated for arbitrary degeneracy and 
without making any approximations in the lepton kinematics: 


A.- 


(6.295 X 10“S"^) (— 

\me 


Ae+ ^ (6.295 X 10-^-^) 




rUf 


E-^ exp(a; - He-lTy + 1 

{x + Q,+ /T,fx{x^-{mjT,fYl^ 

^ exp(x - He+/Te) + 1 


dx; 


dx. 


(43) 

(44) 





29 


Here the Q-values and the energy thresholds are Qe+ ~ 1.293 MeV and = rrie for 
positron capture on neutrons and Qe- ~ —1.293 MeV and = rrin — mp — rUe for electron 
capture on protons. In electromagnetic equilibrium, the sum of the total (kinetic plus rest 
mass) electron and positron chemical potentials is zero, implying that 

Pe+ = -he-- (45) 

We have included neutrino and anti-neutrino captures on heavy nulcei, calculated 
by the prescription given in Fuller & Meyer (1995) and McLaughlin & Fuller (1995). In 
these calculations, we employed Fermi-Dirac distribution functions for the neutrinos and 
antineutrinos normalized by the appropriate average energy and luminosity. The neutrino 
and antineutrino energy luminosities and average energies were taken from the Wilson and 
Mayle output at the relevant time slices. The average neutrino and antineutrino energies 
were £t with straight lines, while the corresponding values for the luminosities were obtained 
using cubic spline interpolation. 

We start our calculations in conditions where weak equilibrium obtains, so that we 
can employ the solution to equation ([Tl| ) as an initial electron fraction. We hold the 
entropy constant throughout our calculation, since we are far from the gain radius when the 
calculation begins. The adiabatic outflow approximation is employed here in keeping with 
our spirit of discerning the basic physics important for Y^. We would expect the entropy of 
an actual outflowing mass element to rise moderately through the regime of weak freeze-out 
(cf. Qian & Woosley 1996). The resulting “non-equilibrium” electron fraction, for three 
different constant entropy trajectories is shown in Figure |^. The resulting electron fraction 
for all the curves follows a downward trend at the earliest times. This trend corresponds 
to the effect of the decreasing importance of the electron and positron capture rates, as 
discussed in section The curve with the highest entropy begins with highest electron 
temperature, since the density history is the same for all of these curves. This also explains 
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the increase in the electron fraction with increasing entropy, since the plasma temperature 
and, hence, the effect of electron and positron capture will be larger at higher entropy. The 
dip in these curves occurs when alpha particles begin to form. This is a smaller version of 
the effect in the “equilibrium” case which was discussed in section ^ It can be seen that 
the smallest “alpha effect” occurs for the case where the electron fraction is closest to 1/2. 

In this graph, heavy nuclei form starting at an epoch between, tpb ~ 0.64 s and 
tpb ^ 0.66 s. After the heavy nuclei form, the system contains mostly heavy nuclei with a 
few percent of the mass fraction in protons. Antineutrinos are capturing on protons, and 
this tends to decrease the electron fraction. However, counteracting this effect are neutrino 
captures on heavy nuclei, which tend to increase the electron fraction. The result is the 
almost flat curves seen at late time in Figure]^. Figure ^further demonstrates the effect of 
neutrino and antineutrino capture on heavy nuclei for a trajectory with entropy s/k k. 40. 
This hgure shows the T/ which is obtained, when captures on heavy nuclei are included 
(upper curve) and not included (lower curve) in the calculation. The captures on heavy 
nuclei represent a relatively small effect, which does not affect the resulting electron fraction 
by more than 1%. 

We hnd that the value of the electron fraction is close to T/q throughout the time 
period where the system is falling out of weak equilibrium. Although, the formation of 
alpha particles, neutrino and antineutrino capture on heavy nuclei, and electron (positron) 
capture on free nucleons can have signihcant leverage on the hnal electron fraction, the 
ratio of the antineutrino to neutrino capture rate on free nucleons at the time of freeze out 
from weak equilibrium has the most influence. However, the relative leverage on W of the 
effects considered here can vary with different outflow conditions. 
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7. Conclusion 

In this paper we have given an in depth treatment and analysis of the evolntion of the 
electron fraction in neutrino-heated supernova outflow, including detailed treatment of the 
effects of nuclear composition changes. Our study has concentrated on the time tpb ^ 1 s. 
The evolution of is this epoch may be quite important for models of the light p-Process 
and the neutron number N ^ 50 r-Process nuclei (Hoffman et ah 1996a; Fuller & Meyer 
1995). 

We have employed hts to the detailed neutrino and antineutrino energy spectra 
from the Mayle and Wilson supernova calculations. We hnd that these detailed spectra 
are necessary for computations of weak rates to the level of accuracy in W which may 
be required for nucleosynthesis considerations. However, we hnd that hydrodynamic 
wave-induced huctuations in the ratios of neutrino and antineutrino spectral parameters 
with time produce signihcant excursions in Yg. We hnd that the rates of electron and 
positron capture on free nucleons can also be important for computing the evolution of the 
electron fraction. The charged current weak rates freeze out from equilibrium at a time 
when the electron and positron capture rates may still have some inhuence on W- During 
this time period, the ehect of these rates is to increase the electron fraction. We have given 
detailed calculations of the “alpha ehect” - the increase in the electron fraction caused by a 
changing alpha particle mass fraction. Our results indicate that the alpha ehect can play a 
very signihcant role in setting W- 

We have employed numerical calculations of nuclear composition changes in nuclear 
statistical equilibrium, coupled with hts to density and velocity of outhow histories 
from the Mayle and Wilson results, to compute the evolution of W- These calculations 
explicitly include all charged current weak interaction processes, including neutrino and 
antineutrino capture on heavy nuclei. We follow the evolution of W through the epoch of 
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weak equilibrium freeze out. These calculations show that the combination of neutrino and 
antineutrino capture on heavy nuclei and antineutrino capture on free protons tend to keep 
Ye constant. The results presented in this paper are meant to illustrate several different 
variables and processes which can alter the electron fraction in post core bonnce snpernova 
ontflow. Clearly, more sophisticated models of nentrino transport and hydrodynamic 
outflow than those employed here would be necessary to actnally predict Ye in a reliable 
fashion. We believe, however, that the effects described here will always play the major role 
in setting Ye. 

We wish to acknowledge usefnl discussions with Y.-Z. Qian, B. S. Meyer, S. E. Woosley, 
and R. D. Hoffman. We also wish to thank R. Mayle for the use of output from the Mayle 
& Wilson supernova code. This work was snpport by NSF grant PHY-9503384 and a NASA 
theory grant at UCSD, while JRW was snpported by NSF grant PHY-9401636. 
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method 

A.e (S-^) 

Ape (s-^) 

Y 

eo 

blackbody spectra 

39.10 

45.86 

0.4602 

cubic spline interpolation 

35.61 

38.96 

0.4775 

pieces of blackbody spectra 

35.54 

38.69 

0.4787 


Table 1: Comparison of effect of (anti)neutrino spectra. 
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Fig. 1.— Comparison of the Mayle & Wilson numerical transport calculation-derived 
neutrino energy spectrum with a zero chemical potential blackbody spectrum. Circles 
indicate the data points taken from the Wilson and Mayle code. The curve through these 
data points has been £t using cubic spline interpolation. The curve without circles shows a 
blackbody spectrum at a temperature of T^, = 3.15 MeV. 

Fig. 2.— The upper curve is the ratio of the neutrino capture rate on free neutrons to 
the antineutrino capture rate on free protons plotted against time. This curve has been 
smoothed in order to average out fluctuations in the neutrino spectra. The lower curve is 
the ratio of the electron capture rate on free protons to the positron capture rate on free 
neutrons plotted against time, in a mass element moving away from the neutron star. 

Fig. 3.— Plot of the quantity Fep against time. The crosses correspond to points calculated 
from neutrino spectral data taken from the Wilson and Mayle supernova code (crosses). The 
same quantity derived similarly from a different run of this code is also shown (circles). The 
line shows the function used in this paper to derive the quantity Wo- 

Fig. 4.— This figure shows the ratio of the positron capture rate on free neutrons to the 
neutrino capture rate on free neutrons against time for an outflowing mass element. 

Fig. 5.— This figure shows the value that the electron fraction of an outflowing mass element 
would take if the system were in weak equilibrium. The lower curve takes into account 
only neutrino and antineutrino captures on free nucleons. The upper curve also takes into 
account electron and positron captures on free nucleons. In calculating these curves it has 
been assumed that the nucleons are free at all times and have not been incorporated into 
nuclei. 

Fig. 6.— This figure shows the value that the electron fraction would take if the system 
were in weak equilibrium. The lower curve takes into account only neutrino and antineutrino 
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captures on free nucleons, as in Figure The upper curve takes into account electron and 
positron captures on free nucleons, as well as the incorporation of free nucleons into alpha 
particles. 

Fig. 7.— This hgure shows the value of the electron fraction derived from a nonequilibrium 
calculation. The entropy was held constant as the mass element flowed away from the 
neutron star. Three different curves are shown, each with a different value for the entropy, 
as indicated. 

Fig. 8.— Results of a nonequilibrium calculation for the electron fraction are shown in 
which neutrino and antineutrino captures on heavy nuclei are included (upper curve) and 
not included (lower curve). The electron fraction is plotted against time post core bounce. 




















